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Abstract 

Let a polyhedral convex set be given by a finite number of linear 
inequalities and consider the problem to project this set onto a sub¬ 
space. This problem, called polyhedral projection problem, is shown 
to be equivalent to multiple objective linear programming. The num¬ 
ber of objectives of the multiple objective linear program is by one 
higher than the dimension of the projected polyhedron. The result im¬ 
plies that an arbitrary vector linear program (with arbitrary polyhedral 
ordering cone) can be solved by solving a multiple objective linear pro¬ 
gram (i.e. a vector linear program with the standard ordering cone) 
with one additional objective space dimension. 

Keywords: vector linear programming, linear vector optimization, 
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1 Problem formulations and solution concepts 

Let k, n,p be positive integers and let two matrices G G H G and 

a vector G be given. We consider the problem of polyhedral projection, 
that is, 

compute Y = {y e M^| G : Gx -|- Hy > h} . (PP) 
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A point (x, y) G M"" xM^’ is said to be feasible for (PP) if it satisfies Gx+Hy > 
h. A direction {x,y) G M"' x (]R^’\{0}) is said to be feasible for (PP) if it 
satisfies Gx + Hy > 0. A pair is said to be feasible for (PP) 

if is a nonempty set of feasible points and is a set of feasible 
directions. 

We use proj : —)■ MP to denote the projection of a set W C M"’+p 

onto its last p components. For a nonempty set B C RP, conv B is the convex 
hull, and conei? :={Ax| A>0, xG convi?} is the convex cone generated 
by this set. We set cone0 := {0}. 


Definition 1. A pair (WP“, is called a solution to ( |PP[ ) if it is feasible, 
AP°‘ and are finite sets, and 


dir 


Y = conv proj AfP°‘ + cone proj X' 

For positive integers n,m,q,r, let A G P G 

h G be given. Consider the vector linear program 

minimize Px s.t. Ax > b, 


( 1 ) 

\ Z G and 


(VLP) 


where minimization is understood with respect to the ordering cone 

G :={yG M'^l Z^y > O} . 


This means that we use the ordering 

w <c y y — w G G 


Z'^w < Z^y. 


( 2 ) 


We assume kerZ^ := G M'^| Z'^y = O} = {0}, which implies that the 
ordering cone G is pointed. Thus, (§ defines a partial ordering, ff Z is 
the qx q unit matrix, (|VLP ) reduces to a multiple objective linear program. 
This special class of (VLP) is denoted by (MOLP). 

A point X G M"' is called feasible for (VLP) if it satisfies the constraint 


Ax > b. A direction x G M”'\{0} is called feasible for (VLP) if Ax > 0. The 


feasible set of (VLP) is denoted by 


5 := {x G M’"| Ax > b} 


A pair (5P°‘, is called feasible for ( VLP| ) if 5P°‘ is a nonempty set of 
feasible points and is a set of feasible directions. The recession cone of 
S is the set O+S* = {x G M”! Ax > 0}. This means that 0+S'\{0} represents 
the set of feasible directions. We refer to a homogeneous problem (associated 


to (VLP)) if the feasible set S is replaced by 0'''5. 
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For a set X C M”-, we write P[X] := {Px\ x G X}. The set 

V ■= p[5] + C = {y G M'?| 3x G M*", Z'^y > Z^Px, Ax > bj 
is called the upper image of (|VLP ). 


Definition 2. A point x € S is said to be a minimizer for (VLP) if there 
is no V G S such that Pv <c Px, Pv / Px, that is, 

xgS, Px + C\{0}. 


A direction x G M"\{0} of S is called a minimizer for (VLP) if the point x is a 


minimizer for the homogeneous problem. This can be expressed equivalently 
as 

xG(0+5)\{0}, Px 0P[O+5] +C\{0}. 

If (5P“, is feasible for ( VLP[ ) and the sets 5^°', are finite; and if 


conv P[SP“] + cone P[S^^^] + C = V 


( 3 ) 


then (5P“,5'^'’') is called a finite infimizer for (VLP) 


A finite infimizer (5P“,5'^‘’’) is called a solution to (VLP) if its two 
components consist of minimizers only. 


This solution concept for (VLP) has been introduced in [1^. It can 


be motivated theoretically by a combination of minimality and inhmum at¬ 
tainment established in [8]. Its relevance for applications has been already 
discussed indirectly in earlier papers, see e.g. mnEiii]. The solver Ben- 
SOLVE [m |T2] uses this concept. 


2 Equivalence between (PP), (MOLP) and (VLP) 


As a hrst result, we show that a solution of ( |PP[ ) can be easily obtained 
from a solution of the following multiple objective linear program 


mm 


^rp ) s.t. Gx + Hy > h, 

6 y j 


( 4 ) 


where e := (1,..., 1)^. Both problems (PP) and (Q have the same feasible 
set but (Q has one additional objective space dimension. 


Theorem 3. Let a polyhedral projection problem ( |PP[ ) he given. If (PP) is 
feasible, then a solution (5P“,5'^'’’) (compare Definition^ of the associated 
multiple objective linear program S exists. Every solution (SP“, of 
(Q is also a solution of (|PP[) (Definition^. 
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Proof. Since (PP) is feasible, the projected polyhedron Y is nonempty. Thus 
it has a hnite representation. This implies that a solution T = 


of (PP) exists. We show that X is also a solution of the associated multiple 


objective linear program 0. X is feasible for a and its two components 
are finite sets. By ([^, we have 


P[S] = 


y 

-e^y 


y€Y'>= convP[XP“] + coneP[5‘^'''] 


which implies ([^ for C = , i.e., T is a finite infimizer. We have 


P[S] C {y G 


^y = o}, 


which implies that all points and directions of S are minimizers. Thus X 
consists of minimizers only. We conclude that a solution of (0 exists. 

Let S = (PP™, be an arbitrary solution of (|^. By (|^ we have 

convP[SP“] + coneP[5'^'’'] + = P[S] + 


where denotes the nonnegative orthant. We show that 


convP[SP“] +coneP[P'^‘’'] = P[S] 


(5) 


( 6 ) 


holds. The inclusion C is obvious by feasibility. Let y G P[S], then e^y = 0. 
By ([^, y G for B := conv P[SP“]+cone There is u G P and 

c G such that y = v+c. Assuming that c G we obtain e^c > 0. 

But B Y P[S] and hence the contradiction 0 = e^y = e^v + e^c > 0. Thus 
^ holds. Omitting the last components of the vectors occurring in ([^, we 
obtain ([^, which completes the proof. □ 


Of course, (MOLP) is a special case of (VLP). In order to obtain equiva¬ 


lence between (VLP), (MOLP) and (PP), it remains to show that a solution 


of (VLP) can be obtained from a solution of (PP). We assign to a given 


(VLP) the polyhedral projection problem 


compute P = {y G M^| G M" : Z'^y > Z^Px, Ax > 6} . (7) 


Obviously, (VLP) is feasible if and only if 0 is feasible. The next re¬ 
sult states that a solution of (VLP), whenever it exists, can be obtained 


from a solution of the associated polyhedral projection problem 0. This 
result is prepared by the following proposition. The main idea is that non- 
minimal points and directions can be omitted in a certain representation of 
a nonempty closed convex set. 
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Proposition 4. Let a nonempty compact set P C of points and a com¬ 
pact set R C M'J\{0} of directions he given and define V := convP + conei?. 
Furthermore, let C C be a nonempty closed convex cone. IfV-\-C C P 
and 

LnC = {0}, (8) 

where L := O+P n —d'^V is the lineality space ofV, then 

V = conv {V\{V + C\ {0})) + cone {R \ (0+P + C\ {0})) + C. 

Proof. The inclusion D is obvious. To show the reverse inclusion, let U be 
a linear subspace complementary to L. Since V and R are compact, and 
F 7 ^ 0, P is a nonempty closed convex set. Then the set ext(P n U) of 
extreme points of V CiU is nonempty, see e.g. [3 Section 2.4]. An element 
V G ext('P n U) admits the representation 

V = ^ +'^Hir\ 

jeJ i&i 

with finite index sets J and I, G V, Xj > 0 for j G J, Ylj&j 
F G R, Hi > 0 for i G I. Every and r* may be decomposed by means of 
= vjj + and r* = + rf with rfi gU and r\ G L. Therefore, 

^ = Z] + Z + Z + Z • 

ieJ is/ j&J i&i 

'-V-" 

ei 

Because v G U, the last two sums vanish. In the resulting representation 

^ = Z + Z FFu 

j&J i&l 

' -V-" 

=0 

the second sum equals to zero because v is an extremal point. For the same 
reason, vfj = v for every j with Xj > 0 follows. Hence, = v + v^j^. 

Assume that gV + C \ {0}. There exist y GV and c G C \ {0} such 
that = y + c. With decompositions y = yr + Uu and c = cl + cu, where 
cjj 0 (because of condition Q and c 7 ^ 0 ), this leads to 

V + v{ = Vj = yu + yi + Cu + Cl, 

whence 

V = yu + cu + yL + CL - . 

' -V-" 

=0 
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Again, the last part being zero results from v being an element of U. Let 
/i G (0,1) be given. As a linear combination of elements of U, yjj + j^cu 
belongs to U. On the other hand, 

yu + -cu = y+-c -yL — cl g V. 

eP GL 

Hence, yu + ^cu € V CiU. By adding a meaningful zero, a representation 
of V, 


V = yu + CU 

= At (^yu + ^cf/^ + (1 - y)yu, 

as a convex combination of different points of "P H t/ is found, which contra¬ 
dicts V being an extremal point. Thus, ^ V + C \ {0}; and altogether 


ext(Pn[/) CH\(p + C\{0}) + L. (9) 

Now an extremal direction r G 0'^{'P DU) is considered. Since V and 

R are compact, V = convH + coneiZ implies 0'''P = coneiZ (see e.g. m 

Corollary 9.1.1]). Thus, a representation 

jeJ jeJ 

' -V-" 

=0 

by G iZ can be obtained. From extremality of r it follows that 

all 7 ^ 0 with fij > 0 coincide with r up to positive scaling. Such an j is 

taken and, without loss of generality, one may assume yj = 1, i.e.: 


We show that 


G (P\ (0+P + C\{0})) U(C + L). (10) 

Indeed, let G O+P + C \ {0}, that is, = y^ + c for y^ G O+P and 
c G C \ {0}. We need to show that ^ C + L follows. If y^ G L, this is 
obvious. Thus, let y^ 0 L. Consider the decompositions y^ = y^ + y^ and 
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c = cu + CL with j/f}, cu G u, y'l, cl & L, y’^j ^ 0 and, by cu 7 ^ 0. We 
obtain 


r + r{=r^=y'lj + y’l + cu + cl- 


Since r G U, 


r = yu + cu + y'l+CL-r^L=yu + cu- 


From L C O+'P the inclusion 0'*'P + L O'^'P and subsequently y^ G 
O'^P — {yi} ^ is deduced. Moreover, cuGC + LCC + 0+P C O’^P. 
Noted that 0+(PnP) = 0+PnP, both directions and cu are in 0+(PnP). 
Therefore the representation of r = y^+cjj as conic combination of elements 
of 0 +(PnP), in which r was supposed to be extremal, proves equality of 
r, y^ and cu up to positive scaling. This implies r G C + L and G C + L. 
From (10) we deduce 


G cone (P\(0+P + C\{0}))+C + L 


and hence 


0+(P n P) C cone (P \ (0+P + C\ {0})) + C + L. (11) 


Any lineality direction I G L can be represented by a conic combination 
of elements of P: I = Ylj&j ■ For any with //fc > 0 this results in 





j&J\{k} ) 


and therefore in G L. If such an is an element of O'^P + C, then there 
exist r G O^P and c G C with = r + c. This implies 



e-o+v 


and by ([^, c = 0 follows. Therefore, G R \ (O+P + C \ {0}), resulting in 

L C cone (P \ (0+P + C \ {0})) . (12) 

Combining the results ( |12[ ) and 0 yields 

L + 0+(P n P) C L + cone (P \ (0+P + C \ {0})) + C + L 
C cone (P\(0+P + C\{0}))+C'. 


7 



Using Q and the decomposition V = L + \P r\U], wq obtain 
V = L + {Vr^U) 

= L + (conv ext(P (1 U) + 0'^ {V H U)) 

C conv (V\{V + C\ {0})) + cone {R \ (O+P + C \ {0})) + C, 


which proves the claim. 


□ 


Theorem 5. Let a vector linear program (VLP) he given. If (VLP) is fea¬ 
sible, a solution of the associated polyhedral projection problem according 
to Definition^ exists. Let X = be a solution of (^. Assume 

that Q is satisfied and set 

:= {x E {x, y) E y ^ p + C\{0}} , 

:= {xEM’^I {x,y)eX^\ y 0 0+P + C\{0}} . 

Then (5P“, 5'*^^“') is a solution of ( VLP ) in the sense of Definition^ Oth¬ 
erwise, if ([^ is violated, (VLP) has no solution. 

Proof. The existence of a solution of the polyhedral projection problem ([^ 
is evident, because a polyhedron has a finite representation. Consider a 
solution (VP“, of ([^. From Proposition]^ we obtain 

V = conv P[5P“] + cone P[S'^^] + C. 


Hence, (5P“, is a finite infimizer for (VLP). It is evident that 5P“ and 
consist of minimizers only. 

Assume now that (|^ is violated. Take y E L n C'\{0}, then for any 
X G S, Px = Px — y + y, where Px — y E P and y E C'\{0}. Thus, x is not 
a minimizer. □ 


The following statement is an immediate consequence of Theorem 


Corollary 6. A solution for (VLP) exists if and only if (VLP) is feasible 
and (® is satisfied. 


In the remainder of this paper we show how a solution of (VLP) can be 
obtained from an irredundant solution of the associated projection problem 
(PP). A solution (XP“, of (PP) is called irredundant if there is no 

solution (I/P™,!/^'’") of (PP) satisfying 

ypoi c J^poi ydir ^ j^dir 


(UP°',Udir) ^ 


dir\ 
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The computation of an irredundant solution of (PP) from an arbitrary so¬ 
lution of (PP) does not depend on the dimension n. It can be realized, for 
instance, by vertex enumeration in 

Theorem 7. Let a vector linear program (VLP) be given. If ( |VLP[ ) is fea¬ 
sible, an irredundant solution of the associated polyhedral projection problem 
0 exists. Let X = he an irredundant solution of 0. Assume 

that ([^ is satisfied and set 

5P“ := {x EM^I {x,y) E X?™} , 

:={xEM"| (x,y)EX^'^ y ^ L + C\{0}] . 

Then (S'P“, is a solution of (|VLPj). Otherwise, if ([^ is violated, 
dVLPl ) has no solution. 

Proof. By Theorem]^ it remains to show that (S'P™, consists of min- 
imizers only. Assume that x for (x, y) E AP“ is not a minimizer. There 
exists z G V and c E C'\{0} such that y = z + c. Let us denote the elements 
of AP°‘ as 

XP“ = {(x\yi),...,(x“-\y“-i),(x,y)}. 

The point 2 ; can be represented by AP“ and d E 0'’"P, which yields 

ol—1 ol—I 

y — c = z = Aj?/* + Xy + d, Ai,..., Aq_i, A > 0, A* -|- A = 1. 

i=l i=\ 


We set V := ^ ^ consider two cases: (i) For A = 1, we have 

u = 0 and hence c = —d, a contradiction to (ii) For A < 1, we obtain 

a■ ■ 1 

y = J2 ''' ^ +cone proj 

i=l 

= conv ((proj AP“) \ {?/}) + cone proj 

This contradicts the assumption that the solution (AP^jA*^^*^) is irredun¬ 
dant. 

Assume now that (x,y) E A'^'^’ with y ^ L + C'\{0} is not a minimizer. 
There exists 2 ; E O+'P and c E C'\{0} such that y = z-\-c. Let us denote the 
elements of A^'^ by 

= [{x\y^),...,{xP-\yP-^),{x,y)}. 
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The direction z can be represented by which yields 

/3-1 

y - C = z = '^Xiy" + Xy, Xi,..., Xy-i, X > 0. 

i=l 

We set V := ^ O+T* and distinguish two cases: (i) Let A > 1. 

Then (1 — X)y — v e —0~^V n C'\{0}. This contradicts ^ since C C 
(ii) Let A < 1. Then 

A. 1 

y = w+d for 7 i ^ ^ > 0, w -fiy'’ and d := ^ G C\{0}. 

2=1 

We have w 0 —0'*’P since otherwise y € L + C'\{0} would follow. For d 
there is a representation 
/9-1 

d = '^yiy" + IJ,y yi,..., yy-i, y > 0. 

2=1 

The condition y > 1 would imply re = (1 — y)y — d £ —O'^V. Thus we have 
/i < 1 and hence 

. I 

y = Y^ ll_^y^ g cone {y\ ..., /“^} = cone (^(proj X'*'’') \ {y}) . 
i=l ^ 

This contradicts the assumption that the solution (XP“,X‘^^'') is irredun- 
dant. □ 


3 Examples and remarks 


It is well known that (VLP) can be solved by considering the multiple ob¬ 
jective linear program 


minimize Z^Px s.t. Ax > b, 


(13) 


compare, c.L, HZ]. The r columns of the matrix Z £ correspond to 

the defining inequalities of C (or equivalently, to the generating vectors of 
the dual cone). The objective space dimension r of (13) can be much larger 


than the objective space dimension q of the initial vector linear program, 
see e.g. m for a sample application and im for the number of inequalities 


required to describe the ordering cone there. In contrast to (13), with our 


approach, the objective space dimension is increased only by one. 

In the following toy-example with q = 2 we illustrate the procedure. 
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Example 8. Let us consider the following instance of (VLP); 


minimize 



s.t. X € S, 


where the feasible set S is defined by 


S := 


X G 



( 14 ) 


(15) 


Minimization is understood with respect to the partial ordering generated by 
the ordering cone 


C := 




= :ZT 


The solution concept for the vector linear program (14) demands 
putation of a representation of the upper image V = P[S] + C, 
depicted in Figure^ First, we express V as an instance of (PP). 


the com- 
which is 
This step 



Figure 1: The feasible set S and the upper image P[S] + C of Example [sj 
A solution is given by the feasible points and the feasible direction x . 

Their respective image-vectors ,y‘^ and y^ generate the upper image. It can 
be seen that x'^ is not part of a solution, as the image-vector y'^ belongs to the 
ordering cone C and is therefore not a minimal direction. 


pushes the ordering cone into the constraint set, compare Q.' 

compute V = {y ^ \ 3x e S, Z'^y < Fx^ , (16) 
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Now we formulate the corresponding instance of (MOLP), with one addi¬ 
tional image space dimension, compare Q; 


minimize 


^0 0 1 0 

0 0 0 1 

,0 0 -1 -1 



s.t. 


^S, 


(17) 


= :P 


where the feasible set S is the same as in (16), 


i.e. 


/ 

( 1 0 \ 


/ 0 0\ 


y A 

' 

/ \ 


1 -1 


0 0 


-1 


< g ^ri+q 


1 1 

X + 

0 0 

A\ 

-1 


\yj 


-1 -3 


-1 2 


0 


K 

\-3 1 y 


V 2 Ij 


\o) 



Now we consider a solution to which consists of := 

with feasible points 



y A 
1 

/™2\ 

-9 X \ 

-1 

-1 

1 

II 

II 

1 

v 1 y 


v-v 


and 


{y^^y^} with the feasible directions 


fx^\ _ 

-1 

w) 

2 


\oJ 


This means y* is M.^-minimal for i = 


fx^\ _ 

0 

w) 

-1 


\^J 


1,..., 4 and 


P[5] + M:; = convP[5P“] + coneP[5^‘’^] + M:;. 

From Theorem @ we deduce that is also a solution for the poly¬ 

hedral projection problem (16). This solution is irredundant, so the x- 
components of the points that is and x"^, are the points of a solution 

to the original vector linear program (14) by Theorem^ It remains to sort 
out those directions whose y-part belongs to the ordering cone C (compare 
Theorem^: This is the case for the direction y^, as Z'^y^ = (5,0)^ ^ 0. 
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Thus the solution for (14) consists of the feasible points x^,x‘^ and the fea¬ 
sible direction x^ (also compare Figure^: 


SPoi ^ 


0 

-1 


^dir ^ 


1 

-1 


It generates the upper image of ([T^ (Figure^ by means of 
V = convP[5P“] + coneP[5^‘'^] + C. 


Finally, let us demonstrate howV can be obtained directly from V. We start 
with an irredundant representation 

V = conv + cone , 

where 



We rule out those y* where the condition = 0 is violated, whence y® is 
cancelled. Then we delete the last component of each vector and obtain 


V = conv {y\y^} + cone {y^,y^] 


with 


y' = ( i)^y" = (_\)’= (o)’ (2 


In the preceeding example one needs one additional objective, wheras the 
“classical” approach (13) does not require any additional objective. Note 


that every step of the procedure presented here is independent of the actual 
values of g ^ 2 and r ^ q. In the following example we consider the case 
r > g + 1. This leads to an advantage in comparison to the classical method 
(13). The cone C C of the vector linear program (18) has 6 extreme 


directions and a solution is obtained from a solution of a corresponding 


the classical approach, see (13), 6 objectives are required 


multiple objective linear program (19) with only 4 objectives. Note that in 


Example 9. Consider the vector linear program 


minimize Px s.t. Bx > a, x > 0 


(18) 
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with ordering cone C = {y G M^| Z'^y > O} and data 

(I 1 1\ /3\ 

12 2 _ 4 

2 2 1’®“ 4 

V2 1 2/ 


P = 



B = 


Z = 


VJ 



4 1 0 0> 
0 0 14 
2 2 2 2 , 


We assign to (18) the multiple objective linear program 

> d, X >0 


minP 

with objective function 


X 

yVy 

P 


s.t. B 


(19) 


y 

-e^y 


and constraints Bx > a, Z'^y > Z^Px, x > 0. Thus, the data of (19) are 


P = 








(1 

1 

1 

0 

0 










1 

2 

2 

0 

0 

0 


4 







2 

2 

1 

0 

0 

0 


4 

/O 

0 

0 

1 

0 

0\ 


2 

1 

2 

0 

0 

0 


4 

0 

0 

0 

0 

1 

0 

D _ 

-6 

-4 

2 

4 

2 

2 


0 

0 

0 

0 

0 

0 

1 

; ^ — 

-6 

-6 

0 

2 

4 

2 

, a = 

0 

VO 

0 

0 

-1 

-1 

-V 


-4 

-2 

2 

4 

0 

2 


0 







-1 

-2 

-1 

1 

0 

2 


0 







-1 

-3 

-2 

0 

1 

2 


0 







V4 

-6 

-2 

0 

4 

2/ 


VO/ 


A solution to (19) consists o/5P“ = {y^,y‘^,y^} with 


y^ = 









0 


0 


0 


1 

1 

to 

II 

II 

2 

-1 


4 

-4 


2 


1 


0 




V2y 


W/ 



(again y® = (x®, 






/o\ 

/o\ 


f = 


0 

0 

0 

-1 

V2/ 


,f = 


0 
0 
2 
2 

V-1/ 


,f = 


0 

1 

-1 

0 

V1 y 


, = 


0 

0 

1 

0 

voy 


, y = 


0 

0 

0 

1 

Voy 
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One can easily check that , y'^, £ C and y^ 0 C. 

the VLP (18) consists o/5P“ = and 5^^“' = 


Thus, a solution to 
|x®|, that is, 


gpoi _ 





^dir 



The upper image V of the MOLP (19) is given by its vertices 








-1 

2 

-2 

1 

-3 

0 

y = 

1 

, y = 

2 

, y = 

4 


v-v 


\-y 


W/ 


and its extreme directions 




^2\ 



-1 

-5 

2 

-6 

0 

2 

, y = 

-1 

, y = 

1 

v-v 




^0/ 



( 0 ^ 



0 

-8 

1 

-9 

0 

0 

, y = 

0 

, y = 

0 

V-1/ 


V-1/ 


VV 


We sort out y^, as e^y^ / 0, and we delete the last component of each vector 
y^,... ,y^. As a result we obtain the vertices 






2/1 = 2 , = 1 1 2/3 = 0 


and the extreme directions 


0 


^_1N 


2/^ = I -1 1 , 2/" = I 2 1 , 2/" = I 0 1 , 2/1 = I 0 1 , 2/“ 

,-l> 


of the upper image V of the VLP (18). 



In the next example we have used the VLP solver bensolve na da 
in order to compute the image of a linear map over a polytope, which can 
be expressed as a polyhedral projection problem. The solver is not able to 
handle ordering cones C = {0}. Therefore a transformation into (MOLP) 
with one additional objective space dimension is required. 
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Example 10. Let W he the 729-dimensional unit hypercube and let P be 
the 3 X 729 matrix whose columns are the 9^ = 729 different ordered arrange¬ 
ments of 3 numbers out of the set {—4, —3, —2, —1,0,1, 2, 3,4}. The aim is 
to compute the polytope P\W\ := {Pw\ w gW}. To this end, we consider 
the multiple objective linear program 


minimize 


Px \ 
—e^Px ) 


s.t. X G W, 


having 4 objectives, 729 variables, and 729 double-sided constraints. The 
upper image intersected with the hyperplane {y G e^y = O} is the set 
P[VE], see Figure^ 



Figure 2: The polytope P\W] of Example [l^ computed by Bensolve p^[T2] 
via MOLP reformulation. The resulting polytope has 43680 vertices and 26186 
facets. The upper image 7^ C of the corresponding MOLP has 43680 
vertices and 26187 facets. The displayed polytope is one of these facets, the 
only one that is bounded. 

We close this article by enumerating some related results which can be 
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found in the literature. It is known [9] that the parametric linear program 

minc^x s.t. Gx + Hy>b (PLP(?/)) 

is equivalent to the problem to project a polyhedral convex set onto a sub¬ 
space. For every parameter vector y we can consider the dual parametric 
linear program 

max(6 — Hy)^u s.t. G'^u = c, u>0, (DPLP(?/)) 

which is closely related to an equivalent characterization of a vector linear 
program (or multiple objective linear program) by the family of all weighted 
sum scalarizations, see e.g. [5]. 

Fiilop’s seminal paper [6] has to be mentioned because a problem similar 
to Q was used there to show that linear bilevel programming is equivalent 
to optimizing a linear objective function over the solution of a multiple 
objective linear program. 

The book [14| provides interesting links between multiple objective linear 
programming and computation of convex polyhedra. 
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